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The buildup process of the probability density inside the quantum well of a double-barrier resonant 
structure is studied by considering the analytic solution of the time dependent Schrodinger equation 
with the initial condition of a cutoff plane wave. For one level systems at resonance condition we 
show that the buildup of the probability density obeys a simple charging up law, [\P (r) /cj>\ — 
1 — e~ T ' T ° , where <j> is the stationary wave function and the transient time constant to is exactly 
two lifetimes. We illustrate that the above formula holds both for symmetrical and asymmetrical 
potential profiles with typical parameters, and even for incidence at different resonance energies. 
Theoretical evidence of a crossover to non-exponential buildup is also discussed. 

PACS: 03.65,73.40.Gk 
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Since the pioneering work of Esaki and Tsu jy, the 
tunneling in one-dimensional semiconductor heterostruc- 
tures has been the subject of intense investigations 
dUHJil- Resonant tunneling in double barrier (DB) sys- 
tems has received special attention both by its technolog- 
ical applications and by the motivation to clarify the new 
interesting transport phenomena. Among the fundamen- 
tal problems that have appeared on the scene, the charge 
buildup in the quantum well region is considered as one 
of the most important processes since it governs the ulti- 
mate speed of resonant tunneling devices H||] . The need 
of a direct and comprehensive dynamical study of this 
phenomenon has been widely recognized How- 
ever, up to now we are lacking of an exact description of 
the buildup process itself. 

In this paper we provide an exact description of the 
buildup process at resonance condition, within the frame- 
work of the shutter model. This is based on a full quan- 
tum dynamical approach, recently developed by Garcfa- 
Calderon and Rubio JtJ, that deals with the solution of 
the time dependent Schrodinger equation for an arbitrary 
potential V{x) (0 < x < L), with an initial condition of 
a cutoff plane wave confined in the half-space x < to 
the left of an absorbing shutter |j| at x = 0. The sudden 
opening of the shutter at t = allows the wave function 
to interact with the potential. As a consequence, they 
found that the transient solution for the internal region 
may be written as the stationary solution modulated by 
a time varying Moshinsky function plus an infinite sum of 
transient resonance terms associated with the S— matrix 
poles of the problem. 

For the case of a reflecting shutter 
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one can proceed along the lines similar to that discussed 
in Ref. and obtain the solution for the internal region, 

V(x,k;t) = 4>{x,k)M{Q,k;t) - (p*(x,k)M(0, -k;t) 
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where (/>(:r, k) is the stationary wave function and the fac- 
tors T n — 2ku n (0)u n (x) I (k 2 — k 2 ) are given in terms of 
the resonant eigenfunctions u n (x). The index n runs over 
the complex poles k n distributed in the third and fourth 
quadrants in the complex fc-plane. The Moshinsky func- 
tions , as it is well known, are defined in terms of the 
complex error function w(z) [|l(|: M{y q ) = M(0,q]t) — 
w{iy q )/2, where the argument y q is given by 
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and q stands either for ±k or k± n . 

The time evolution of wavefunction ^(x, k; t) in the in- 
ternal region may be described by the expression given 
by Eq. (|J), which involves the contribution of the full 
resonant spectrum of the system. For structures with 
typical parameters, and incidence energies E = h 2 k 2 /2m 
near a resonance energy e„, Garcfa-Calderon and Rubio 
|Q showed that the single resonance approximation for 
the wavefunction ^(x,k;t) is valid from a few tenths of 
the corresponding lifetime onwards. This is the case for 
the present study, since we are not considering the regime 
of very short times (t < r„ = h/T n ), which may require 
the contribution of far away resonances. The single res- 
onance approximation to Eq. (H) is 
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where we have used the fact that the poles in the third 
quadrant fe_ n are related to those of the fourth, k n , by 
k- n = — fc*. This is the one-level expression of the time 
dependent wave function for the description of the dy- 
namics in the internal region. 

In order to exemplify the building up of the probability 
density in the quantum well for incidence at different res- 
onance energies and for different potential profiles, let us 
consider the following two numerical examples. The first 
case corresponds to the symmetrical DB structure with 
parameters: barrier heights Vq = 0.5 eV, barrier widths 
bo = 30 A and well width ujq = 100 A. The resonance pa- 
rameters for the first three resonant states are: e\ = 37.8 
meV, Ti = 0.12 meV; e 2 = 149.2 meV, T 2 = 1.40 meV; 
e 3 = 325.7 meV, T 3 = 8.60 meV. We show in Fig. 1 the 
time evolution of the probability density calculated by 
Eq. (|J) for incidence at resonance, E = e n and fixed po- 
sition x (we have considered values of x near the maxima 
of \<j)(x,k) | 2 as the most natural choice) for the cases: 
n = 1, x = 80 A (solid line); n — 2, x — 48 A (dashed 
line); and n — 3, x — 80 A (dotted line). The second 
example consists of an asymmetrical DB structure with 
parameters: barrier heights Vi = V2 = 0.3 eV, barrier 
widths bi = 30 A and 62 = 100 A, and well width luq = 50 
A. The resonance parameters for the first resonant state 
are Si = 89.1 meV and Ti = 2.4 meV. The time evolu- 
tion of 1^ (a;, k; t) | 2 for incidence at E = £1 and x — 55 
A is also depicted in Fig. [I] (dashed-dotted line). For 
all cases, the probability density \^(x,k;t) | 2 grows up 
monotonically towards its asymptotic value. We can see 
that both the level off and the rate of increase of the 
curves are quite different. However, a common feature, 
not evident in Fig. [I], can be appreciated if we replot 
the normalized probability density (x,k;r) /<p(x,k)\ 2 
as a function of the new variable r, which is now the 
time given in lifetime units (i replaced by Th/T n for each 
curve). As a result, all four curves become indistinguish- 
able among them as depicted in Fig. ||. We can see that 
the full establishment of the stationary situation is pre- 
ceded by a transient in which the probability density is 
built up inside the quantum well with a unique character- 
istic curve. Thus, there must also exist a characteristic 
transient time constant To that governs the buildup pro- 
cess, with the same value (in lifetimes) for all cases. The 
observed regularity and the fact that it holds for both 
symmetrical and asymmetrical cases and at different res- 
onances, is a manifestation that the buildup process in 
one-level systems is governed by a simple law. In what 
follows we shall be concerned to find the analytic expres- 
sion of the actual buildup law and the exact value of the 
characteristic transient time To. 

For sharp (R n = e n /T n ^> 1) and isolated reso- 
nances we know that the stationary wavefunction <f>(x, k) 
can be written as the one-term expression |TT[ <f){x, k) = 
2iku n (0) u n (x) /(k 2 — fc 2 ), thus the factors iT n and iT_„ 
appearing in Eq. (0) can be identified as cp(x, k) and 



—cj)*{x,k) respectively. By expressing formula (^) in life- 
time units, it can be seen that y q depends only on the 
ratio R n = s n /T n , and not on the particular values of 
the resonance parameters e n and T n . For q — ±fc„ and 
q = ±k n , y q reads: y ±kn = Te^ 4 [(Rn - i/2) r] 1/2 , 
and y±fe» = =pe - " r / 4 [(i?„ + i/2) t] 1 ^ 2 , respectively. For 

the cases q — ±k we have y±k = T e ~ I7r ^ 4 [RnT] 1 ^ 2 (since 
E = e n ). From the above considerations and the well 
known symmetry relation [Q, M{y q ) = e y i — M(—y q ), 
applied to the Moshinsky functions M(y^) and M(yk n ), 
we obtain a convenient representation for the probability 
density, 

I* (x, k; r)| 2 = \4>{x, k)\ 2 (1 - e-^ 2 ) 2 + A(r), (5) 

where A(r) stands for the remaining terms, which involve 
the square modulus of the Moshinsky functions and sev- 
eral interference terms. It is not difficult to convince 
oneself that, for very large times, expression @ pos- 
sesses the correct asymptotic behavior, i. e. as r — » 00, 
1^ (x, k; t) \ 2 goes into the stationary probability density 
\<fi(x, k)\ 2 . This follows directly from the presence of the 
decreasing exponential e~ T / 2 , and from the fact that each 
of the Moshinsky functions involved in A(r) can be repre- 
sented by a series expansion consisting of inverse powers 
of r. In Ref. it was shown that M(y q ) has the asymp- 
totic expansion M(y q ) = a\/y q + a 2 /y 2 + a 3 /y* + ■ ■ ■ 
for large values of the variable y q provided that — n/2 < 
arg(y q ) < tt/2. By inspection of our expressions for y q , 
we can see that the above inequality holds for the three 
cases involved in A(r) (q = —k, —k n , and —A;*), and as 
a consequence A(t) — » when r — * 00, as expected. It 
is important not only that both e _r / 2 and A(r) go to 
zero, but also the fact that their rates of decrease are 
quite different, leading to important consequences on the 
nature of the buildup. In particular it reveals that there 
exist exponential and non-exponential contributions to 
the buildup mechanism, as we shall see later. In view 
of the series expansions considered above, A(r) also con- 
tains inverse powers of y q , that is, inverse powers of the 
product (-R^r) 1 / 2 . This means that A(r) can be vanish- 
ingly small even for r equal to a few lifetimes, provided 
that Rn 3> 1. In fact, there exists a finite time interval in 
which A(t)/|0| 2 is negligible compared to e~ r / 2 , leading 
to the following exponential buildup law: 

|*(r)/0| = l- e -/ 2 . (6) 

This simple formula reproduces successfully the pre- 
dicted values of expression (^). For comparison, we have 
included in Fig. || a plot of the normalized probability 
density calculated from (^|). The corresponding curve is 
indistinguishable from all the other curves, showing in 
particular that A(r) has an exceedingly small contribu- 
tion in the relevant time interval for all of our numerical 
examples. We see that formula (pi) does not depend ex- 
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plicitly on the potential profile parameters nor the reso- 
nant state, this explains why all the numerical examples 
illustrated in Fig. || share the same curve, despite the 
fact that they correspond to different situations. Note 
that in the exponential regime the buildup law becomes 
identical to the charging up law of a capacitor in an 
RC— circuit: Q{t)/Qq = 1 — e~ r / rc , where Q is the 
asymptotic charge, and tq = RC is the capacitive time 
constant. This is relevant, because we find in the lit- 
erature capacitor-like models used to describe quantum 
tunneling properties in DB structures, such as the charge 
buildup and its implications on the speed limit on reso- 
nant tunneling devices [^|Jl^1. According to Eq. (H), the 
transient time constant tq of our "quantum capacitor" 
is always two lifetimes, and is a characteristic feature of 
one-level systems. On the experimental side, it is worth 
to mention that measured values of escape times of the 
order of 2K/T n has been reported by Sakaki et al. ||, 
arguing that, at coherence conditions, "the buildup time 
and the tunneling escape time are roughly the same" [Q. 
An important remark is that the condition R n 3> 1 is not 
so restrictive since it is satisfied for most of the resonant 
structures with typical parameters. In fact we carried out 
a systematical study (not shown here) and found that for 
values of R n from 10 onwards this condition is satisfied. 

In order to show the existence of deviations from the 
exponential regime, we shall examine the contributions 
arising from A(r). The explicit calculation of A(r) in 
terms of y~k, D-k n , and y-k* may result a too involved 
task, however, for the purpose of our discussion, it is suf- 
ficient to realize that the dominant term is proportional 
to an oscillatory function of r modulated by the factor 
t -1 / 2 . We know that at very long times the exponential 
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term goes to zero faster than r 
Therefore, there must exist a critical time r onse t at which 
e~ T / 2 and A(t)/\4>\ 2 are comparable. Such a critical 
time defines a crossover from the exponential to a non- 
exponential regime of the buildup process. In the exam- 
ples depicted in Fig. ^ the non-exponential contributions 
to (^(t) /<p\ 2 are overwhelmed and cannot be appreci- 
ated due to the scale of the graph. However, if we plot 
the logarithm of the difference 6(t) = |1 — |W (r) /cf>\\ ver- 
sus r, [using Eq. (Q)] the transition from the exponential 
to the non-exponential regime is clearly appreciated, see 
Fig. [| In this figure, the exponential regime can be 
identified by the straight line with slope —1/2, extend- 
ing over a few lifetimes until it reaches the onset of the 
nonexponential buildup, T onset , which depends on R n . 

It is interesting to note the similarity of the results 
depicted in Fig. || to the behavior of the survival prob- 
ability found in studies of the phenomenon of quantum 
decay Jl^Jl^l , which also exhibits this transition with an 
oscillatory structure at such crossover. We believe that 
this striking resemblance is not a simple coincidence but 
rather a manifestation of the existence of a more pro- 
found link between both phenomena. In fact, the sur- 



vival probability may also be expressed in terms of the 
Moshinsky functions Jl3| , which are, in our expressions, 
the key ingredients for the time evolution. These find- 
ings open up new questions about the common features 
in both processes, for example those about the existence 
of deviations from the exponential buildup also at early 



times, as it occurs in the decay process 15 1. Such anal- 
ysis requires the contribution of far away resonances to 
the transient solution, and is deferred to future work Jl6| . 

Summarizing: (i) We have accomplished the first an- 
alytic derivation of the actual buildup law in resonant 
tunneling structures. It was based on general proper- 
ties of the solution of the Schrodinger equation, without 
any assumptions on the potential profile, except that it 
is finite and support well defined resonances, (ii) We 
have shown the existence of both exponential and non- 
exponential contributions to the buildup process. (Hi) 
The exponential regime is characterized by a transient 
time constant whose value is exactly two lifetimes, (iv) 
We have illustrated that formula (||) describes very ac- 
curately the exponential buildup for a great variety of 
situations: it works very well for different potential pro- 
files, and is valid not only for the "ground state" (n = 1), 
but also for "excited states" (n > 1). 

We thank G. Garcfa-Calderon for useful discussions. 
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FIG. 1. This graph illustrates the building up of the prob- 
ability density for the two examples discussed in the text. 
In the symmetrical system we show the time evolution for 

(x, k; t) | 2 for incidence at: si = 37.8 meV (solid line), 
e 2 = 149.2 meV (dashed line) and e 3 = 325.7 meV (dotted 
line). The fixed positions x are: 80 A, 48 A and 80 A, respec- 
tively. In the asymmetric case e\ = 89.1 meV (dashed-dotted 
line) at x = 55 A 

FIG. 2. Evolution of |* (r) /<j>\ 2 as a function of r for the 
same cases of Fig.l. All four curves become identical. For 
comparison, a plot of |* (r) /<f>\ 2 using Eq.(6) is also included. 
The resulting values are also superimposed on the other curves 
illustrating that they obey a simple exponential law with a 
unique transient time constant of two lifetimes. 

FIG. 3. Exponential and non-exponential contributions to 
the buildup. We plot the logarithm of the difference of 
S(t) = |1 — |\I>(t)/0|| versus r, for the states n = 1 and 
n = 3 of the symmetrical system considered in the text. The 
corresponding values of the ratios R n = e n /V n are show in 
the figure. The linear behavior corresponds to the exponen- 
tial regime, and the deviations from it appear after a certain 
transition time r onset . 
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